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Abstract. Our goal in this paper is to extend the theory of quasi-exactly 
solvable Schrodinger operators beyond the Lie-algebraic class. Let V n be the 
space of ra-th degree polynomials in one variable. We first analyze exceptional 
polynomial subspaces M C T n , which are those proper subspaces of Vn invari- 
ant under second order differential operators which do not preserve Vn- We 
characterize the only possible exceptional subspaces of codimension one and 
we describe the space of second order differential operators that leave these 
subspaces invariant. We then use equivalence under changes of variable and 
gauge transformations to achieve a complete classification of these new, non- 
Lie algebraic Schrodinger operators. As an example, we discuss a finite gap 
elliptic potential which does not belong to the Treibich-Verdier class. 

1. Introduction 

Our main focus of study in this paper is the class of second order linear ordinary 
differential operators T with analytic coefficients which preserve a finite dimensional 
vector space 

M =span{/i,...,/„}, 
where the fi(z) are polynomials in some coordinate z. For a given M., we define 

V 2 {M) = {Te V 2 \TM C M} 

to be the vector space of second order differential operators that preserve M.. 

The class of operators which belong to the space T>2(M) for some M. is of 
considerable interest from both a mathematical and a physical perspective since 
it encompasses, after a suitable change of variable, almost all 1 the known exam- 
ples of potentials for which the Schrodinger equation is either exactly solvable, or 
quasi-exactly solvable in the sense that part of the spectrum together with the cor- 
responding eigenfunctions can be determined algebraically, ft is also the basis of 
much of the modern theory of orthogonal polynomials. 

Historically, the Lie algebra b^K) has played a rather important role in the 
development of the subject, since the vector space V n of univariate polynomials 
of degree at most n is an irreducible s[2(M)-module for every positive integer n. 
This implies that a second-order differential operator which can be expressed as 
a polynomial in the generators of s^R) acting on Vn will automatically preserve 
V n . Since the 5I2 action on V n is irreducible, Burnside's Theorem ensures that the 
converse is also true, i.e. every operator that preserves V n belongs to the enveloping 
algebra of 5(2 . 



^We should point out that the hydrogen atom potential, and more generally, the class of 
Natanzon exactly-solvable potentials do not fall directly into this framework, but that they can 
be recovered in terms of invariant spaces of polynomials through a slightly different approach. 
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When the operator is put in Schrodinger form by a suitable choice of independent 
variable and gauge, the resulting potentials, which are referred to as Lie algebraic 
potentials, have been widely studied in the literature. They include the class of 
Poschl- Teller, Morse, sextic anharmonic, and Lame potentials. One of the lim- 
itations of the Lie algebraic approach is that there is no a-priori reason for it to 
encompass all the exactly or quasi-exactly solvable potentials. Furthermore, it does 
not provide a systematic way to detect the presence of multiple algebraic sectors 
in the spectrum, corresponding to different subspaces which are preserved by the 
same operator. 

The general point of our paper is that the class of (quasi-)exactly solvable poten- 
tials is considerably richer than the Lie algebraic class. We are thus lead to introduce 
the important concept of an exceptional subspace M. , defined as a proper subspace 
of V n generated by n-th degree polynomials such that T>2{M) </- T>2(V n ). Excep- 
tional subspaces lead naturally to exceptional potentials which are not Lie- algebraic. 
In this paper we characterize all exceptional subspaces of co-dimension one, denoted 
as Xi subspaces (Theorem l4.ll) , and we provide explicit basis for T>2 (Xi ) (see Propo- 
sition^^!). Most remarkably, both the new Xi operators and the Lie algebraic 
ones admit an a-priori characterization at the level of the potential (Theorems 13.51 
and I4.11|) , expressed as a linear or a quadratic constraint over the residues of the 
quotient q(z)/p(z) of the coefficients of the operator T = p(z)D zz + q(z)D z + r(z). 
One of the most useful aspects of this a-priori characterization is that it quickly 
allows to decide which potentials are non-singular, and which ones have multiple 
algebraic sectors. Making use of this characterization, and the projective action of 
SL(2, R) on V n , we classify all the new Xi potentials into canonical forms, providing 
in each case the potential form in the physical variable with the constraints over 
the coefficients and the gauge factor. 

As an example of the new potential forms, in Section we analyze in greater 
detail a non-singular elliptic potential on the line which is a deformation of the 
Lame potential, admitting four algebraic sectors. This potential is a finite gap 
potential which does not belong to the Treibich-Verdier class [19]. The existence of 
these types of finite-gap potentials was proved in [20] but no results were given on 
their solvability. We prove that the potential is quasi-exactly solvable and we show 
how to compute the algebraic eigenfunctions that correspond to the band edges 
(similar to the Lame polynomials). These generalized finite-gap potentials have 
been studied recently in [17]. 

Summarizing, our work illustrates once again the fact that the representation 
theory of s^R) is not needed in order to construct any of the quasi-exactly po- 
tentials, even the Lie algebraic ones, and that reliance on the purely Lie algebraic 
setting imposes unnecessary limitations to the theory. 

Our paper is organized as follows. In Sectional we recall some basic facts on the 
transformation properties of second-order ordinary differential operators, and the 
equivalence problem to operators in Schrodinger form. Section [3] is devoted to the 
Lie algebraic class of Schrodinger operators and potentials which arise from finite- 
dimensional irreducible representations of sfe (K) by first order differential operators. 
The classification of Lie-algebraic potentials reproduces the classifications in [10, 
21] with the added feature that potentials with multiple algebraic sectors are also 
classified. Section 0] studies exceptional subspaces of co-dimension one, and the 
second order differential operators that leave them invariant. A characterization of 
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the new operators at the level of their potential invariant is also given. The full 
classification of Xi potentials into canonical forms under the SL2 projective action 
is given in Section[3] An Xi elliptic finite gap potential is treated in Sectional with 
an explicit calculation of its eigenfunctions. Finally Section sums up the results 
and outlines future developments. 




2. Second order operators. 

In this section we consider the general equivalence problem for second order 
differential operators in one variable. 

Let z be a coordinate on a 1-dimensional domain, and let D z be the unique first 
order operator such that D z [z] — 1. Consider a general second order operator 

T=p(z)D zz + q(z)D z +r(z), (I) 

where p, q, r are analytic functions of z. We will call p(z), q{z) and r(z) the compo- 
nents of T relative to the coordinate z. We first observe that the components of T 
are determined by knowing the action of T on three linearly independent functions. 

Proposition 2.1. Let fk{z), k — 1,2,3 be three smooth, linearly independent 
functions, and set gk{z) = T[fk(z)]. The components of a second order operator T 
are the unique solutions of the following linear equation: 

<K fl ft 
& & h 
\f'{ f'z 1z 

The components of a second order differential operator undergo two very natural 
transformations: changes of coordinates and gauge transformations. 

A smooth, invertible change of coordinate z = C( w ) transforms the components 
according to 

KCH ) 
C'H 2 

_ g(cM) p(gw))c"(w) 

q[ ' C'H C'H 3 ' 

f(w) = r(CH), 

where p, q, f denote the components relative to w. 

Wc define a gauge transformation to be the conjugation of an operator by a 
multiplication (zero-order) operator: 

f = e a{z) Te- a{z} . (3) 

The components of the transformed T are given by 

p(z)=p{z), (4) 

q{z) = q(z)-2p(z)a'(z), 

f(z) = r(z) - q{z)a'{z) + p(z){a'(z) 2 - a"(z)) 
We will say that a second order operator S is in self-adjoint form if 

S = p(z)D zz + ^p'(z)D z + v(z). (5) 
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Every second order operator T can be transformed by a gauge transformation into 
self-adjoint form S. Given T with components p(z),q(z) and r(z) one must take 
o{z) as 

i r 



aiz) = 2j W) dz - 



v(z) = r(z) - -q\z) + -p"(z) - (q(z) - \p\z)) (q(z) - §,/(*)) 



In this way the gauge transformation Q transforms T into self-adjoint form with 

= ^)-kw + V'W-i^# + i^. (0) 

2 4 4 p(z) 16 p(z) 

Since © is invariant with respect to gauge transformations (0J we shall refer to 
v(z) as the potential invariant or simply, the potential of operator T. 
We will say that an operator H is in Schrodinger form if 

H = -D xx + u(x) 

If p(z) < on the domain of interest, then an operator S in self-adjoint form can be 
transformed into Schrodinger form by a change of coordinate. The required change 
of coordinate is the inverse of the function x = £(z), given by 

«*)= ' 



The potential u(x) is then given by 

u(x) = v(((x)). 

Summarizing, an arbitrary second order differential operator T can be trans- 
formed into a Schrodinger operator by a change of variable and a gauge transfor- 
mation. We note that this result is no longer valid for differential operators in two 
or more variables, where the equivalence problem is much harder [14]. This is the 
main obstacle to the extension of quasi-exact solvability to multi- variable operators. 

3. The Lie algebraic class. 

As mentioned in the Introduction, our main purpose is to study differential opera- 
tors that preserve finite dimensional polynomial subspaces. More specifically, given 
a polynomial subspace M. we want to give a basis of T>2(A4), the set of second order 
differential operators that leave M. invariant. We then use the transformations de- 
scribed in the previous Section to construct and classify Schrodinger operators with 
finite dimensional invariant subspaces of functions. When the invariant polynomial 
space is 

V n (z) = (l,z,z 2 ,...,z n ) (7) 

the resulting class of potentials is the well known Lie algebraic class of quasi- 
exactly solvable potentials, [10,21]. In this Section we reproduce the classification 
of Lie algebraic potentials using an approach which is different from the usual one 
based on the generators of SI2 realized as first order differential operators. Our 
approach allows to characterize which operators are quasi-exactly solvable at the 
level of the potential invariant of the self-adjoint form. The advantage of this 
characterization is that a certain symmetry in the condition on the parameters 
becomes apparent, which allows the same Schrodinger operator to be algebraized 
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in more than one way, thus leading to multiple algebraic sectors. Examples of this 
phenomenon had already been discussed in the literature [4,6,9, 13], but in this 
Section we provide a complete classification of the Lie-algebraic potentials that 
admit multiple algebraizations. The reason to revisit the well known classification 
of Lie-algebraic potentials is that the same approach can be extended to other 
polynomial modules, which we deal with in Section 

3.1. Second order operators that preserve V n - The structure group for the 
equivalence problem for differential operators is the infinite-dimensional group gen- 
erated by diffeomorphisms and gauge transformations. If an operator preserves 
a finite dimensional vector space of polynomials then the structure group of the 
equivalence problem is reduced to SL2. 

The n + 1 dimensional vector space of polynomials 

V n (z) = (l,z,z 2 ,...,z n ) 

can be made into an irreducible SL2 representation. The irreducible action of SL2 
on V n is given by 

f(w) = ( 1W + 5) n f(aw)), f(z)eV n (z), (8) 

where 

JW + 

is a fractional linear transformation corresponding to an element of SL2 . The above 
action is an irreducible multiplier representation [15] of SL2, which is isomorphic to 
the unique n + 1 dimensional irreducible representation of SL2. The corresponding 
infinitesimal generators of the SI2 action are the following first order operators 

T_ = D z , T = zD z - |, T+ = z 2 D z - nz, (10) 

which leave invariant the n + 1 dimensional space V n {z). 

The above SL2 action extends naturally to an action on second order differen- 
tial operators. Let T and T be second order differential operators related by the 
following gauge transformation: 

f[f(w)} = ( 1 w + 5) n T[f(CM)], (11) 

where f(z) and f(w) are related by iJSJ, and z = ((w) is the fractional linear 
transformation shown in ©. Applying J2J and the transformation law for the 
operator components is seen to be 

p(w) = ( 1 w + 8) 4 p(CH), (12) 
q(w) = ( 7 ui + 5) 2 q(((w)) - 2(n - 1)7(7^ + <5) 3 P(CH), 
f(w) — r(C(ui)) — 717(71/7 + 5)q(£(w)) 

+ n(n - l)7 2 ( 7 u; + 5) 2 p(((w)), 

while the transformation law for the corresponding potential invariants © is simply 



v(w) — v(C,(w)). 



(13) 
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The inversion transformation £(u>) = —1/w plays a significant role in the subsequent 
analysis. We will therefore say that a second order operator T is the inversion of 



operator T if the two operators are related by 

f[f(w)} = w n T[f(-l/w)}, f(w) = w n f(-l/w). (14) 
Specializing H12(l to the case of an inversion gives the following transformation rules: 

p(w) = w A p(—l/w), (15) 

q(w) = w 2 q{-l/w) - 2(n - l)w 3 p(-l/w), (16) 

f(w) = r(—l/w) — nwq(—l/w) + n{n — l)w 2 p(— 1/w) (17) 
As well, a calculation shows that 

p'(w) = w 2 p'{-l/w) +4w 3 p(-l/w). (18) 



Since acts irreducibly on V n , Burnsidc's Theorem [21] ensures that a second 
order operator T preserves V n if and only if it is a quadratic element of the en- 
veloping algebra of the SI2 operators shown in (|10[l . Thus, the most general second 
order differential operator T that preserves V n can be written as 

T= c iJ T i T J+ b * T * + c *> ( 19 ) 

i,j=±,0 i=±,0 

where = Cji, bi and c* are constants 2 . For this reason, an operator that preserves 
V n {z) is often referred to as a Lie-algebraic operator. Note that Burnside's Theorem 
does not apply to other subspaces M. generated by polynomials, and therefore an 
operator T <G T>2{M) is not necessarily Lie algebraic. 

In the rest of the Section we present a different characterization of operators 
that preserve V n . Let 

n 

denote the infinite-dimensional vector space of all polynomials. 

Proposition 3.1. A second order operator preserves V{z) if and only if all its 
components are polynomials in z. 

Proposition 3.2. If a second order operator preserves V n [z) for some n>2, then 
it preserves Viz). 

Proof. Suppose that T = p(z)D zz + q(z)D z + r(z) preserves P n (z). We apply 
Proposition 12 . II with l,z,z 2 and obtain 

r(z)=T[l], q (z)=T[z]-zT{l], p(z) = \ T[z 2 } - zT[z] + ^T[l]. 

□ 

Proposition 3.3. Let T be a second order differential operator and T its inversion. 
The operator T preserves V n (z), n > 2 if and only if T preserves V(z) and T 
preserves V(w). 



In the Lie algebraic approach it is more difficult to determine the dimension of T>2 {Vn ) than 
using the so called direct approach [7]. The family 1191 has a priori ten free parameters, but there 
is one relation among them coming from the quadratic Casimir of sb, thus d\mT>2( r P n ) = 9. 
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Proof. If T preserves V n (z), then it preserves V{z) by Proposition 13.21 By SL2 
covariance, as expressed by the transformation law 114(1 . the transformed operator 
T preserves V n {w) and hence V(w). 

Conversely, suppose that T and T preserve V(z) and V(w), respectively. Let 
f{z) e V n {z) be given and set 

g(z) = T[f(z)}, f(w) = w n f(-l/w) e V n (w), g(w) =f[f(w)]. 

By assumption, g(z) € V(z) while g(w) € V(w). However, because of covariance 

g(w) = w n g(-l/w), 

which implies that degg(z) < n. Therefore, T preserves V n (z). □ 

We can now give a basis for T>2(V n ). There are other, more efficient ways to 
arrive at this result, but we pursue this particular approach because it will allow 
us to study T>2(M) where M. is a more general polynomial space. 

Proposition 3.4. If a second order operator T = p(z)D zz + q(z)D z +r(z) preserves 
V n {z), then necessarily degp{z) < 4. Indeed, T preserves V n (z) if and only if it is 
a linear combination of the following 9 operators: 



Proof. Using ((15(1 - 1(17(1 a calculation shows that the inversions of the operators 
shown in ((20(1 - 1(23(1 have polynomial coefficients. Therefore, by Proposition the 
above 9 operators preserve V n (z). 

Let us now prove that the above operators span T>2(V n ). Let T = p(z)D zz + 
q(z)D z + r(z) be a second order operator that preserves V n (z). By Proposition ^. 31 
the components of the inversion of Tare polynomials in z. An examination of 1(15(1 
then shows that, necessarily, degp(z) < 4. Hence, subtracting a linear combination 
of the operators l(2"n jl -l|22 )l . we obtain a first order operator T\ — q\{z)D z + ri(z) 
that preserves V n (2). An examination of ((16(1 shows that, necessarily, degqi(z) < 2. 
Hence, by subtracting a linear combination of the operators shown in 1(23(1 . we arrive 
at a multiplication (zero-order) operator that preserves V n (z). By fI7jl. such an 
operator can only be a constant. □ 

3.2. The characterization of Lie-algebraic potentials. Following [22] we re- 
call that the Lie-algebraic character of an operator can be manifest at the level of 
its potential invariant. We shall say that v(z) is a Lie-algebraic potential if there 
exists a second order Lie-algebraic operator T — p(z)D zz + q{z)D z + r(z) such that 
v(z) is the corresponding potential invariant given by ©• 

Below we recall the classification [10,21,22] of Lic-algebraic potentials based on a 
case-by-case analysis of the various configurations of the roots of the fourth degree 
polynomial p(z). We begin by describing the generic case: 

Theorem 3.5. Let p(z) be a J^th degree polynomial with distinct roots Zi, and let 
k\,k2, kz, &4 be constants satisfying 



z 4 D zz - 2(n - l)z 3 D z + n(n 
z 3 D zz - 2(n - l)z 2 D z + n(n 

z 2 D z — nz, zD z , D z ,\ 




(20) 
(21) 
(22) 
(23) 



^1 + ^2 + ^3 + ^4 = — n - 1. 



(24) 
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Then, the rational function 

v{z) = -f2(k 2 -±)^ + \. (25) 

i 6 Zi 

is a Lie-algebraic potential. Conversely, let T — p(z)D zz +q(z)D z + r(z) be a second 
order operator such that p{z) is a J^th degree polynomial with distinct roots Zi; let 
v(z) denote the corresponding potential, as given by and set 

If T preserves V n (z), then equations (|24|) and l|25|l hold. 

Proof. Suppose that p(z) satisfies the stated assumption, that k\ , ki , kz , ki satisfy 
(|24|l . and that v(z) has the form shown in (|25J) . Without loss of generality, p(Q) ^ 0. 
Set T = p(z)D zz + q{z)D z + r(z), where q(z) is the polynomial defined by the 
equation 



and where 



q(z) ^ 2t, + l 



Ku^ 1 ///.^ , Uq-P'(z)) 2 1 P'{zf 



(27) 



r(z) = ±q'(z) - -p"(z) 



2 w 4 w 4 16 p(z) 

-£(tf-&) f^ + A- 

2—1 

In this way, v(z) is related to the just-constructed operator T by formula (JSJ). Taking 
the limit as z — > Zi in (|27|) gives (|2(j|l . A calculation shows that the residues of r(z) at 
z — Zi all vanish, and hence r(z) is a polynomial. Therefore, by Proposition 13. II T 
preserves V(z). Next, let T = p(w)D ww +q(w)D w +f(w) denote the inversion of T, 
and let v(w) = v(—l/w) be the corresponding potential. Making the substitutions 
z = —1/w, Zi = —1/wi, and using equation (|18f) we obtain 

0H = -^( fei 2 -i) -^Ua , (28) 



where 

A 

1 S P'(Wi) 



Ao = A-5> 2 -^)^. (29) 



2—1 



Applying lfl3|) and l(TB|l gives 



gH 1 g(-!/w) , s 1 , . 

~ r T^ = ~2 ? i/ N +(2-2n) — . (30) 
j»(w) «r p(— l/u>) w 



QUASI-EXACT SOLVABILITY IN A GENERAL POLYNOMIAL SETTING 



Equations JSJ and {23 then give 

f4=f2-2n-^ + l)U + E ? T ±1 < 31 > 

P[W) \ L — ' / If Z — ' W — Wi 

V ' \ 1=1 / 2—1 

z— 1 

Hence, mutatis mutandi, the above argument can be applied to conclude that T 
preserves V(w), and by Proposition 13. 31 the operator T preserves V n {z). 

Conversely, suppose that T = p(z)D zz + q(z)D z + r(z) preserves V n (z). Let v(z) 
be defined by 10 and fcj be the constants defined by l|26|l . We denote by 

»/ s /x 1 i, 1 // / x 1(q — p'(z)) 2 1 p'(^) 2 

= r « - + - i^y^ + Te^i (33) 

1 

the difference of expressions (0 and (|25[) . By Proposition 13.21 and r(z) are 
polynomials. A calculation shows that A(z) has vanishing residues at z = z t , and 
is therefore a polynomial. Let T — p(w)D ww + q(w)D w + f(w) be the inversion 
of T, and v(w) = v(—l/w) its potential invariant. Let A(w) denote the difference 
of v(w) and the expression in the right-hand side of (|28H . Repeating the above 
argument, we conclude that A(w) is a polynomial in w. Making the substitutions 
z = —l/w,Zi = —1/wi in and using l(T%)l we obtain that A(w) = A(— 1/z), 
which is only possible if both sides are constant. Hence, v(z) has the form, up to a 
constant term, shown in (|25|) . Since both q(z) and q(w) are polynomials, equation 
(I30|l implies that degq < 3. Hence, equation l|2bj) implies equation l|27|l and since 
q(w) is a polynomial, then equation (|31() implies i|24ll . □ 

3.3. The generalized Lame potentials. Suppose that p(z) has four distinct 
roots. In describing the generic Lie algebraic potential, no generality is lost if we 
assume that one of the roots is at infinity, i.e., that degp = 3. One can further 
specialize the form of the above potentials by stipulating 

1 

zi = 1 , Z2 = 0, Z3 = 1, 2:4 = 00, < m < 1. 

m 

Specializing Theorem 13. 51 we obtain the following. 

Proposition 3.6. Let S be an operator in self-adjoint form and suppose that 
p(z) = 4z(l — z)(mz — to + 1). The corresponding potential v(z) is Lie-algebraic if 
and only if 

(m — 1) 

777Z — 777+1 



v(z) = 4 (h 2 - ±) 777 z - (fc 2 2 _ " -±l - + (34) 



such that 

fci + fc 2 + &3 + &4 = — n — 1. (35) 



10 



DAVID GOMEZ-ULLATE, NIKY KAMRAN, AND ROBERT MILSON 



3.4. Multiple algebraic sectors. It has long been known that some exactly solv- 
able potentials possess multiple algebraic sectors. For example, the harmonic oscil- 
lator possess one algebraic sector corresponding to even eigenfunctions, and another 
corresponding to odd eigenfunctions. The existence of multiple algebraizations for 
the Lame potential has been noted in [1,4] 

We will say that a second-order operator T possesses multiple SL2 algebraizations 
if there exists a local coordinate z such that T preserves both the vector space 
V n {z) and the vector space <\>{z)Vn{z). A distinct, second algebraic sector arises 
if 4>(z) is not a polynomial, leading to more algebraic eigenfunctions. In some 
cases, when 4>( z ) is a polynomial, it happens that §{z)Vh{z) C "P n ( z ) an d n0 new 
algebraic eigenfunctions arise. We say then that the operator T possesses nested 
SL2 algebraizations. 

Multiple SL2 algebraizations have been discussed in [1,4,9], but a full classifi- 
cation of potentials admitting multiple SL2 algebraizations was not known. Here 
we provide the full classification thanks to the explicit potential form shown in 
equations l)24 [l -l)25 [l . In the generic case, when p{z) has simple roots, multiple SL2 
algebraizations arise from the parameter symmetry 

4 

k = (-l) Pi ki, n = n+^(l-(-l)f%, Pl e{0,l}. 

z=l 

in the formula for v(z) shown in l|25|) . The possible parameter symmetries are given 
by the following 

Proposition 3.7. Let n be an integer. The above h is an integer if and only if 2ki 
is an integer for some i, or 2(fcj ± kj) is an integer for some i 7^ j . 

Consequently, there are 4 cases in which multiple algebraic sectors arise. 

(M2a) If ki £ \ Z for exactly one i, there are two algebraic sectors. 

(M2b) If there is no i such that fcj G 4 Z, but ki±kj £ ^ Z for exactly two choices 

of i =/= j, then there are two algebraic sectors. 
(M4a) If ki £ \ Z for exactly two values of i, there are four algebraic sectors. 
(M4b) If ki £ j Z for all i, there are four algebraic sectors. 
(M8) If h £ \ 7L for all i, there are eight algebraic sectors. 

We observe that if the root structure of p(z) is degenerate, i.e. if degp(z) < 2 
the opportunity for multiple algebraizations is lessened, but some possibilities still 
remain. The full classification of multiple SL2 algebraizations is given in Tables IT151 
As above, the classification is achieved by examining all the possible symmetries in 
the parameters of the potential. 

Note that here we restrict ourselves to multiple SL2 algebraizations, in which the 
invariant polynomial spaces are of the same type. However, other multiple alge- 
braizations exist [8] in which the invariant polynomial spaces are of a different type, 
namely V n and an exceptional polynomial space, like those described in Section 0] 

3.5. Classification of Lie-algebraic potentials. We finally describe the classi- 
fication of Lie-algebraic potentials. The novelty with respect to previous classifica- 
tions is that the potentials which admit multiple algebraic factors have also been 
classified. 
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We will refer 3 to the potentials in proposition ^ . 61 as type I . lip(z) has a multiple 
root, no generality will be lost if we assume that this root is located at infinity. 
Normalizing the value of the other root(s) to various convenient values, and taking 
limits of (|34[1 and (|35|) . as appropriate, we obtain the remaining cases. 

3.5.1. Type I potentials. 

p(z) = 4z(l — z)(mz — to + 1), z — cn(x|m) 2 , to € (0, 1), 

v(z) = 4 (k^ - i) to z - (fc 2 2 - i) ^— (36) 

mz — to + 1 

, /. 2 1 x (^-1) , /, 2 is 1 

+ ( fc 3 - TsJ + ( fc 4 - Te) 7^Y' 

n = — (fci + &2 + fa + ki) — 1, 

fi(x) = dn(ai|m) C2 cn(x|m) C3 sn(x|TO) C4 , a = 2k t + 1 /2. 

Taking fc 3 ,/c4 = ±1/4, makes the potential non-singular on the interval < z < 
1. The resulting class is referred to as the generalized or the associated Lame 
potentials [9]. The ordinary Lame potentials are obtained by setting k 2 = ±1/4. 
The possibilities for multiple algebraizations are indicated in Table ^ Requiring 
that the potential be non-singular eliminates many of the possibilities. 



4> n — n Conditions 



_l + i)- 2fc 2^-2fc 3 ( z _ 1 )-2'=4 


2(fc 2 


+ ki + k 4 ) 


fell > \k 2 + k 3 +k 4 


z~ 2k3 (z- l)- 2ki 


2(k 3 


+ k 4 ) 


ki + k 2 \ > |fc 3 + k 4 




2(fei 


+ k 3 + ki) 


k 2 \ > \ki + k 3 + k 4 


\ m ' 


2{k 2 


+ k 3 ) 


k\ + ki\ > \k 2 + k 3 




2(fei 


+ k 2 + k 3 ) 


k 4 \ > \ki +k 2 +k 3 


(z-l + ±)- 2k Hz-l)- 2k * 


2(fe 2 


+ k 4 ) 


ki + k 3 \ > |/c 2 + ki 




2(fei 


+ k 2 + k 4 ) 


k 3 \ > \ki + k 2 + ki 


( z - 1 + J_)-2fc2 


2k 2 




k!+k 3 + ki\ > \k 2 




2(fei 


+ k 2 ) 


k 3 + k 4 \ > \ki + k 2 


z -2k 2 


2fc 3 




fel +fe + M > |/B3 




2(fci 


+ k 3 ) 


fc 2 + > |fcl + fc 3 


(z - l)~ 2ki 


2fc 4 




fci + fc 2 + fc 3 | > \k 4 




2(fei 


+ k 4 ) 


k 2 +k 3 \> \ki + ki 


1 


2fei 




k 2 +k 3 + ki\ > |fci 



Table 1 . Multiple algebraizations for type I potentials 



It is worth mentioning that the class of type I Lie algebraic potentials described 
above coincides with the well known class of Treibich-Verdier elliptic finite gap 
potentials, [13,18,19]. 



We borrow the nomenclature for the root structure of 4th degree polynomial from the Penrose- 
Petrov classification of the 4-dimensional Weyl tensor. 
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3.5.2. Type II potentials. 

p(z) = 4z(l — z), z = — sinh 2 x, 

v(z) = Ah 2 z 2 - 44(4 + 24) * - fc32 ~^ ~ ki [lf ' ( 3? ) 
n = -(£ 2 + fc 3 + fc4)-l, 

2 1 

fi(x) = cxp(— l\ sinh a;)(sinhx) C3 (coshx) C4 , Cj = 2fc; + — . 

There are multiple algebraic sectors if one of the three parameters 4,^3,^4 is a 
half-integer. If all three parameters are half-integers then there are four algebraic 
sectors (possibly nested). 

The non-singular potentials arise when k 3 — ±1/4, and in this case the sign of 
&3 governs the parity of the eigenfunctions. Thus, there is one type of non-singular 
type II potentials with multiple algebraic sectors, (line 1 of Table [5J . 



4> 


h — 


n 


Conditions 


z- 2k *(z~\y 2k ± 


2(k 3 i 


-fc 4 ) 


141 > |fc 3 4-fc 4 | 


z -2k a e ~2£ lZ 


2(4 4 


h) 


N > |4 + fc 3 | 


(z - l)- 2k * e~ 2eiz 


2(4 4 




N > \h + k 4 \ 




24 




\k 3 + k 4 \ > |4| 


z -2k 3 


2k 3 




|4 + fc 4 | > \k 3 \ 


( z _ l)-a** 


2ki 




\£ 2 + k 3 \ > \ki\ 



Table 2. Multiple algebraizations for type II potentials 



3.5.3. Type D potentials. 

p(z) = — z 2 1 z = exp(a;), 

v(z) = £ l 2 z 2 - 244 2 + £ 3 2 z- 2 + 244 z~ x , (38) 

n = -{£ 2 +£ 4 )-l, 

fj,(x) = exp(4e x - £ 3 e~ x +44 

These potentials are never singular. There are two algebraic sectors whenever 
4,4 e |Z. 



<p h — n 


Conditions 


e (2e 3 )/z z -2i, 2fi 


141 > 141 


e 2 ^ z 2£ 2 


141 > |4| 



Table 3. Multiple algebraizations for type D potentials 
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3.5.4. Type Z potentials. 

p(z) = -{z 2 + 1), z = sinh:r (39) 

,22 oee 244* + 4 2 -4 2 + i/4 

v(z) = £i z z z -244 z H '—, 

z z + 1 

rc=-(4+4)-l, 

jit(a;) = exp {4 sinhx + 4 tan _1 (sinha;)} (coshx)^ 4 . 

These potentials are never singular. There are 2 algebraic sectors if 4,4 are half- 
integers. 



4> n — n 


Conditions 


e- 2 ^ z 24 


141 > |4| 


e -2£ 3 ta„- 1 ( z )( z 2 + 1 )-£ 4 21a 


141 > |4| 



Table 4. Multiple algebraizations for type Z potentials 



3.5.5. Type III potentials. 

p(z) = -4z, z = x 2 , (40) 

v(z) = £ 3 2 z 3 - 444 z 2 + 4 (4 2 + 44)* + (fc 4 2 - to) z-\ 
n = -(4 + fc 4 ) - 1, 

^(x) = exp{-±4a; 4 + £ 2 x 2 } x Ci , c 4 = 2fc 4 + i. 

These are the well-known QES sextics. They are non-singular if and only if fc 4 = 
±1/4, in which case the sign of fc 4 governs the parity of the cigenfunctions. There 
are two algebraic sectors if and only if 4,&4 & \ Z. Thus, a non-singular QES 
sextic possesses only one Lie-algebraic sector. 

4> h — n Conditions 

e ^ 3 z 2 ~2i 2 z 2h | fc4 | > 

z- 2k * 2fc 4 |4| > N 
Table 5. Multiple algebraizations for type III potentials 



3.5.6. Type N potentials. 

p(z) = -l, z = x, (41) 

v(z) = (l lZ + l 2 ) 2 z 2 + {l x z + 4X4* + 4), 

4 = -2 - 2n, 

fj,{x) = exp {\£\x z - \iix 2 - \l 3 x] 
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To obtain normalizable eigenfunctions, one must take i\ = 0, £2 > 0. The result is 
a harmonic oscillator potential. 

4. Exceptional subspaces. 

Let Ai C V n be a finite dimensional subspace generated by polynomials. We 
define M. to be an exceptional subspace ofV n if T>2{Ai) <f_ T>2(V n )- The significance 
of exceptional subspaces is clearly that they are the only ones that lead to quasi- 
exactly solvable potentials which do not belong to the Lie-algebraic class. For 
brevity we will denote by Xfc an exceptional subspace of co-dimension k. In the rest 
of the paper we study and classify all QES potentials whose invariant subspace is 
Xi. 

The analysis [9, 16] of polynomial subspaces spanned by monomials brought to 
light two special subspaces: 

£„(z) = (l,z 2 ,...,z"), (42) 

£ n (z) = (l,z,z 2 ,...,z n - 2 ,z n ). (43) 

These two subspaces £ n and £ n are SL2-equivalent, since 

4H = w n £ n (-l/w) = {w n f(-l/w) : f(z) e £ n (z)}. 

The following result shows that £ n and £„ are not only special within the class of 
spaces generated by monomials, but that they are essentially the only Xi subspaces. 

Theorem 4.1. Let A4 C V n be an n- dimensional polynomial subspace (not nec- 
essarily spanned by monomials), and suppose that M. 7^ "P n — !■ Then, A4 is an 
exceptional subspace if and only if M is SL2- equivalent to £ n . 

This theorem limits the only interesting codimcnsion-one polynomial subspaces for 
which new QES potentials can be constructed. The proof is rather lengthy and 
shall be given elsewhere. 

We devote the rest of the paper to the classification of the new Xi potentials. 
For this purpose, we need to give a basis of T>2(£ n ) and then describe the potentials 
that correspond via © to such operators. First, we inquire about the effect of an 
SL2 transformation on £ n . The effect of an inversion transformation, z = —1/w, 
was described above. More generally, let us set 

£(z) = (l,z 2 ,z 3 ,...), (44) 

£^ h (z) = (a(z -b)-l,(z- b)\ (z - bf, ...), (45) 

£ a n \z) - (a(z -b)-l,(z- b) 2 , (z - b) 3 , ...,(z- b) n ) , (46) 

4(z) = (l,z,z 2 ,...,z"- 2 ,z"-az"- 1 ), 

and note that £^ b = £ a b n V n . 

Proposition 4.2. The infinite- dimensional vector space £ a - b consists of polynomi- 
als f G V that satisfy by the following first order constraint: 

f'(b)+af(b) = Q (47) 

Proposition 4.3. The subspaces £^ b and £% are Xi. To be more precise, let 
z = C(w) = {aw + 13)/ (~fw + S) , aS — f3j = 1 
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be a fractional linear transformation. If a ^ 0, we have 

(TW + S) n S n (C(w))=£Z\w), (48) 
where a = —nwy, b = — (3/a. If a = 0, then 

( 1 w + 5r£ n (((w))=£«(w), (49) 

where a = —nS/'f. 

Proof. In the first case, using aS — f3j = 1, we have 

(jw + 5) n £ n (((w)) = <( 7 u> + 5) n - k {aw + f3) k : k = 0, 2, . . . , n) 

= ((-a(w -b)+ n) n ~ k {w - b) k : k = 0, 2, . . . ,n> , 

= £ a n\w). 

The last equality is true because (— a(w — b) +n) n ~ k (w — b) k satisfies <|47[1 for k = 
and for k > 2. If a — 0, then 

( 7 w + <T) n £ n (CH) = ((«;- a/n) n - k : k = 0, 2, . . . , n) 
= (1, ti;, w 2 , . . . , w"" 2 , / - aw™" 1 ) 

Since £% ,b and £^ are SL2 equivalent to £ n , by Theorem 14. II thev are Xi. □ 

We describe explicitly the relationship between exceptional subspaces related by 
an inversion w = —1/z: 

Proposition 4.4. For b 7^ 0, the Xi subspaces are related by an inversion 

u n £Z> b {-l/u)=£Z'\u), where a = b(n + ab), 6 = -1/6. (50) 
For b — 0, we have 

u n £%°(-l/u)=£- a (u). (51) 

Proof. Let z = {aw+ (3) / (p/w+S), aS—(3j = 1 be a fractional linear transformation. 
Suppose that a, (3 7^ 0. With to = —1/u, we have z = (fiu — a)/(5u — 7). Setting 

b = —(3/a, a = —nor], b — aj ' (3, a — —n(3S, 

we have, by Proposition 14.31 

( 7 w + 5) n £ n (z) = £%< b (w), 
{5u- 1 ) n £ n {z)=£i\u). 

From this, lf5U|) follows. 

To prove (|51J) . suppose that z, w, u are as above, but that (3 — 0. Since nj/8 = 
na^j = —a, we have 

{ 1W + 5) n £ n {z) = £^{w), 
(6u- 7 ) n £ n (z)=£- a (u). 

□ 
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Our main interest in this Section is to study second order operators T that 
preserve £^ ,b (z) but do not preserve V n (z), which we shall denote as Xi operators. 
These are the only non-Lie algebraic operators (up to SL2 equivalence) that preserve 
a codimension one polynomial subspace, and thus lead to new QES potentials. 

The characterization of T>2(V n (z)) was straightforward: an operator that pre- 
serves V n {z) must map polynomials to polynomials and the same must be true 
for all SL2-related operators. However, the characterization of T>2{£n' b (z)) is more 
complicated: the operators in question must preserve the first order condition i|47|) 
in an SL2-covariant manner. This is addressed by Propositions 14.51 and |4~§1 which 
should be regarded as the Xi analogues of Propositions 13. II and 13.31 respectively. 

Proposition 4.5. A second order T = p(z)D zz + q(z)D z + r(z) preserves £ a ' b (z) 
if and only ifp(z) is a polynomial, and there exist polynomials qo(z) and ro(z) such 
that 

q(z) = q (z)+p>(z)-^, (52) 

z — 

r(z) = r (z)- 2 -^±, (53) 

9o(6) = 2ap(6), (54) 

r' (b)=aq' (b) + ap"(b)-a 2 p'(b). (55) 

The above Proposition follows directly from the following two Lemmas. 

Lemma 4.6. A second order operator T preserves £ a,b {z) if and only if it is a 
linear combination of 

D zz + [2a- D z - (56) 
\ z — J z — b 

{z-b)D zz + (a(z-b)-\)D z , (57) 

(z-b) k D„,k>2, (58) 

(z-b)(D z +a), (z-b) k D z , k>2, (59) 

1, (z-b) k , k>2. (60) 

Proof. Showing that the above operators preserve £ a ' b (z) is a a straightforward 
calculation. We prove the converse. Since £ a ' b (z) — £ a '°(z — b), no generality will 
be lost by restricting to the case 6 = 0; the general case follows by applying a trans- 
lation transformation. Suppose that T is a second order operator that preserves 
£ a <°(z), so that T[az - 1], T[z\ T[z 3 }, T[z 4 ] are all polynomials. By Proposition 
12. II and a calculation, it follows that T is a linear combination of 

z~ x D zz , z- 2 D Zl z~ x D z , z~\ z- 2 

and an operator with polynomial components. Hence, T = T% + Ti, where T2 is a 
linear combination of the operators (|56|l - (|60[) . and T\ is an operator of the form 

Ti = D zz + [q H 1 T D z + nz H 1 T . 

z \ z z z / z z z 

However, T\ applied to az — 1, z 2 , z 3 , z A must yield polynomials that satisfy (|47ll . 
which is only possible if T\ = 0. □ 

Lemma 4.7. A second order operator satisfies conditions 1)52(1 - Q55J1 if and only 
if it is a linear combination of the operators 1)56(1 - (|60[1 . 
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Proof. One direction of the claim is a straightforward calculation. Regarding the 
converse, let T be an operator that satisfies (|52|l - (|55|l . Subtracting the appropriate 
linear combination of i|56|) - i|58|) from T yields a first order operator that also satisfies 
(I52|l - (|55() . so it suffices to prove the claim for a first order T , i.e., we can assume 
p(z) = 0. In this case, conditions (|52|l and (|53[1 assert that q(z) and r(z) are 
polynomials, while conditions (|54|) Q55JI reduce to 

<?(&) = 0, r'(6) = ag'(6). 

By inspection, the first order T in question must be a linear combination of the 
operators shown in (|59|I - H6Q(I . as was to be shown. □ 

Proposition 4.8. // a second order T preserves £%' b (z), n > A, then T preserves 
the full £ a ' b {z). 

Proof. The proof of Lemma |4.6I applies, mutatis mutandi, to show that if T pre- 
serves £ n' h (z) for n > 4, then T is a linear combination of the operators (|56(l - (|60|l . 
and these operators preserve £ a,b (z). □ 

Proposition 4.9. Let n > 4 and a,b be given. Suppose that 6^0 and set a = 
b(n + ab), b = —l/b. A second order operator T preserves £^' b (z) if and only ifT 
preserves £ a ' b (z) and its inversion preserves £ a,b (w). 

Proof. Suppose that T preserves £^' b (z). By Proposition 14.41 the inversion of T, 
call it T, preserves £^' b (w). The desired conclusion now follows by Proposition ^. 81 
To prove the converse, let f(z) G £n' b (z) be given and set 

g(z)=T[f(z)}, f(w)^ W n f(-l/w)e£^' b (w), g(w) = f[f(w)}. 
By assumption, g(z) S £ a,b (z) while g(w) £ £ a,b {w), and by covariance, 

g (w) = w n g(~l/w). 

This implies that degg < n, and therefore, T preserves £^ b (z). □ 

Proposition 4.10. // a second order operator T preserves £^ b (z), then necessarily 
degp(z) < 4. Furthermore, a basis of T>2{£n b (z)) is given by the following seven 
operators: 



Jl 


= (z- 


bfD zz 


- 2(n - l)(z - b) 3 D z + n(n - l)(z - bf , 


(61) 


Jl 


= {z- 


bfD zz 


-{n-l){z-bfD z , 


(62) 


Js 


= {z- 


b) 2 D zz 




(63) 


Ji 


= {z- 


b)D zz - 


V{a{z-b)-l)D z , 


(64) 


J 5 


= D ZZ 


+ 2 (a 


f \ 2a 

l) D » 1 ' 

z — b J z — 


(65) 


Je 


= (z- 


b) (z- 


b--)D z -n(z-b), 


(66) 


Jj 


= I. 




(67) 



Proof. There are several ways to obtain the above basis. The proof below is based 
on SL 2 covariance, c.f., the proof of Proposition l3.4l See [9] for a different approach. 

Without loss of generality, we assume that b 0. The above operators satisfy 
conditions (|H2J)-l(E3), and hence preserve £ a ' b (z). Using lTr^ - l(T7|l . a calculation 
shows that the inversions of these operators obey an d hence preserve 
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£ a,b (w), where a, b are given by (|50|) . Therefore, by Proposition 14.91 the above 7 
operators preserve £%' b (z). 

Let us now prove that the above operators span T>^{£^ ' (z)). Let T = p(z)D zz + 
q(z)D z + r(z) be a second order operator that preserves £^ b {z) and let T = 
p{w)D ww + q(w)D w + f(w) denote the inversion of T. The latter preserves £% ,b (w) 
by Proposition 14. 41 By Propositions 14. 51 and 14. 81 both p(z) and p(w) are polynomi- 
als, whence by i[15[l. necessarily, degp < 4. Hence, by subtracting an appropriate 
linear combination of the operators J\, J2, J3, J4, J5 we obtain a first order operator 
T\ = q\{z)D z +r\{z) that preserves £® ,b (z). Let T\ = q\(w)D w +r\(w) denote the 
inversion of T\ , where by (|16[1 (|1 T|) we have 

qi(w) — w 2 qi(—l/w), fi(w) = r(—l/w) — nwq(—l/w). 

By Propositions 14.51 and l4~8l qi(z),qi(w),ri(z) and ri(w) are polynomials; hence 
deggi < 2 while degri < 1. By Lemma |4.7I T\ must be a linear combination of 
(z — b)(D z + a), [z — b) 2 D z , and a constant. The requirement that fi(w) be a 
polynomial forces 2\ to be a multiple of Jq plus a constant. □ 

We shall refer to the potential invariant v(z) of an Xi operator as an Xi potential. 
We can now state the main theoretical result of the paper: the characterization of 
Xi potentials. It turns out that it is easier to describe Xi potentials than Xi 
operators as evinced by the following 4 

Theorem 4.11. Let p{z) be a fourth degree polynomial with distinct roots Zi, and 
let b be a constant such p{b) 7^ 0. Let kx,k%, k% and k 4 be constants such that 



k\ + fc 2 + fc 3 + k 4 = — n, (68) 
= 0. (69) 



Then, the rational function 



is an Xi potential. The corresponding Xi operator preserves £n ,b (z), where 

(71) 



Eki 
b~ 

1=1 



Conversely, let T — p{z)D zz + q(z)D z + r(z) be a second-order operator such that 
p{z) is a fourth degree polynomial with distinct roots Zi- Let v(z) denote the corre- 
sponding potential, as given by (©, and let 

2p'(zi) 2 

If T is an Xi operator preserving £^ b (z), then p(b) 7^ 0, and equations 1)68(1 ■ <|69[1 . 
f70}l. (|7T|) all hold. 



4 It is instructive to compare the characterization of Xi potentials given in Theorem l4,lll with 
the corresponding characterization of Lie algebraic potentials given in Theorem 13^51 
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The proof of the theorem proceeds via a number of Lemmas. The following Lemma 
identifies those operators that are both Lie-algebraic and preserve £%' b (z). Above, 
we adopted the convention that an operator cannot be both Xi and Lie-algebraic. 
Thus, some Lie-algebraic operators preserve an Xi subspace, but by definition, an 
Xi operator must not preserve V n (z). 

Lemma 4.12. Let T = p(z)D zz + q(z)D z + r(z) be a second order operator such 
thatp(z) is a fourth degree polynomial with distinct roots z^. IfT preserves £°" b (z), 
thenp(b) = if and only ifT also preserves V(z). Similarly, ifT preserves £%' b (z), 
then p(b) = if and only ifT also preserves V n {z). 

Proof. Suppose that T preserves £ a,b (z). If T also preserves V{z), then by Proposi- 
tion l3.ll q(z) is a polynomial. Hence, by equation (|52|l of Proposition ^. 51 we must 
have p(b) = 0. Conversely, suppose that p(b) = 0. By 1(52(1 and 1)53(1 both q(z) and 
r{z) are polynomials. Hence, T preserves V{z), by Proposition 13.21 

Next, suppose that T preserves the finite-dimensional £ ^ ,b (z), and let T denote 
the inversion of T. If T preserves V n {z), then the above argument shows that 
p(b) — 0. Conversely, suppose that p(b) — 0. A repeat of the above argument shows 
that T preserves V{z). A calculation based on (|15l) - (|16fl shows that p(b) = where 
b = —1/6. By Proposition 14.41 T preserves £^ ,b (z), and therefore the transformed 
components q(w) and r(w) are polynomials in w. Proposition l3.3l thcn implies that 
T preserves V n (z). □ 



Lemma 4.13. Let p{z) be a fourth degree polynomial with distinct roots Zi and 
qo(z) be a polynomial. Let a and b be constants such that equation 1(54(1 holds and 
p(b) 7^ 0. Let T — p{z)D zz + q(z)D z + r(z) be a second order operator where q(z) is 
given by l(52|) and where r(z) is a rational function. Let ki be defined by ((72(1 and 
let A(z) denote the difference of the right-hand side expressions in © and 1(70(1 . 
Then, T preserves £ a ' b (z) if and only if A(z) is a polynomial such that 



Proof. We begin by listing some consequences of the above assumptions. Define 
r (z) so that equation ifSHjl holds. Substituting ifSl^ into the given expression 
for A(z) gives 




(73) 



A{z)=r (z) 




2ap(z) 1 
z-b ~ 2 



l (q (z)(z-b)-2p(z)) 2 1 p\zf 
4 p{z)(z-bf 16 p(z) 




) 



(z - 6) 2 




4 



+ 



-A. 
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Rearranging, and using 1(54(1 . the above expression can be written as 

A(z) = r (z) - U(z) - \p"{z) - 2a ( E^lzM) (74 ) 



2^' y ' 4" y ' V z-b 
<*>(*) -qa(b) iP '(z)-p'(b) 2 fp{z)- iAh)-(:-h ]f ,' ( U) 



z — b z — b \ (z — b) 

1 P'{zf lq a (z) 2 > J2( k 2 l_j P'(zi) A 



16 p(z) 4 p(z) 
Applying the identity 

f p\zf 



z=b i=l V ' 



x P( 2 ) y 

using ((5411 , and canceling gives 

A'(6) = (/ (6) - ap"(b) + a 2 p'(b) - aq' (b)) - £ (76) 

Note that ((55(1 describes the vanishing of the bracketed term in the right hand side 
of the above equation. As well, (f72|) and (f5*2|l give 

k - qo(Zl) (77) 

We also note that, by equation l|74|) . A(z) has vanishing residues at z — %{. 

Now, suppose that T preserves £ a,b (z). By Proposition ^. 51 ro(z) is a polynomial 
that satisfies equation (J5HJ. Therefore, by equations itTH) (tTBTl (t77|) . A(z) is a 
polynomial such that condition ((73(1 holds. Conversely, suppose that A(z) is a 
polynomial such that 173(1 holds. Since A(z) has no poles, equations 174(1 177(1 
imply that ro{z) is a polynomial. Finally, equations 173(1 and 1(76(1 imply that 155(1 
holds. Therefore, by Proposition ^. 51 T preserves £ a ' (z). □ 

The following Lemmas establish additional transformation rules for the inversion 
z = — l/io. 

Lemma 4.14. Let p{z) be a fourth degree polynomial with distinct roots Zi =/= 0; 
let b be a constant such that p(b) ^ 0; let v(z) be defined by l(7U(l : set z = —1/w, 
b = —1/6, Wi — —1/zi, v(w) = v(—l/w). Then, 

< w ) = — ( TT 2 1 -2^{ k i ~ Iff) — - + A °' ( 78 ) 

(w — by w — b ~~T w — w l 



whe 



As well, 



y ktp>{ Zi ) ^y kjp'i^) (?9) 
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Proof. Applying the above substitutions and 1|18|) yields 

2p(b) p'(b) _ 2p(b) p'(b) | 2p(b) p'(b) 
(z~b) 2 z-b (w-b) 2 w-b b 2 b ' 

p'(zi) = p'(Wj) p'(Wj) 

z — Zi w — Wi w% ' 
from which equations l|78|) (|79|) follow immediately. □ 

Lemma 4.15. Let p(z) be a fourth degree polynomial with distinct roots z% 7^ 0; let 
p(w) be the polynomial defined by (|15|1 and let q(z) and q(w) be rational functions 
related by Ijl6|l : let qo(z) be defined by equation l)52[l . and let qo(w) be defined by 
the analogous equation: 

q( w )^q ( w )+p'( w )-^H, (80) 
w — b 

where w = — 1/z, b = —1/6. Then, (|72[1 holds and qo(z), q~o(w) are both polynomials 
if and only if (|68f) holds and if qo{z) obeys 



Qo(z) 



2p(z) ^ z - 

Proof. Using equations l|15|) (|16|l l|18|l l|52|) . a calculation shows that 

Qoiw) 1 q (-l/w) 2n 



(82) 



p(w) w 2 p(—l/w) w 

Suppose now that l|72|l holds and that qo(z) and qo(w) are both polynomials. By 
equation this is only possible if degqo < 3. Equation (|72|l gives equation (I77|l . 
which gives equation (|81|l . 

Conversely, suppose that equations (j68(l and (|81|l hold. Taking the residues of 
the latter at z = Zi gives equation 177(1 , which gives equation (|72|l . Using equations 
(|68[) l|81|) (|82|l and making the substitutions z = —l/w,Zi = —1/wi gives 



A k I A \ 

\ ^ / . \ , \ -1 



go(w) = ^ 

2p(w) ^ w — u). 



= E^^- ( 83 ) 

•f— ' W — Wi 

Therefore, qo(w) is also a polynomial. □ 

Proof of theorem \4-.ll\ Suppose that p(z) and b satisfy the stated assumptions, 
that k\, &2, fc 3 , /C4 are constants satisfying equations I|tj8fl and i|t)9|l . and that u(z) 
has the form shown in (7UJ. We assume, without loss of generality, that p(0) 7^ 0. 
Define the polynomial qo(z) by (|81|l and set q(z) so that l|52|l holds. Define the 
constant a by equation (|71|l : evaluating l|81|) at z — b then shows that equation 
(|54|l holds. Consider the operator T = p(z)D zz + q(z)D z + r(z), where r(z) is 
defined so that © holds. By construction, A(z) = 0, where A(z) denotes the 
difference between the right-hand sides of © and (|70J) . Hence, by Lemma f4. 131 T 
preserves £ a - b {z). Let T = p(w)D ww + q{w)D w + f(w) be the inversion of T, and 
let v(w) = v(—l/w) denote the corresponding potential. By Lemma 14.141 v(w) has 
the form shown in 17U|) with w,p(w) 1 Wi = —l/zi, and 6 = —1/6 replacing z,p{z), Zi, 
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and b , respectively. Define qo(w) so that equation l|8Q|l holds. By Lemma [4.151 
qo(w) is a polynomial having the form shown in Q8H|I . Hence, mutatis mutandi, we 
can apply the above argument to conclude that T preserves £ a ' b (w), where 

g = £ ,3- = (-h- - ki) = b(n + ab). (84) 

Therefore, by Proposition 14.91 the operator T preserves £„ ,b (z). By Lemma [4.121 
T does not preserve V n (z), and is, therefore, an Xi operator. 

Conversely, suppose that T = p(z)D zz + q z D z + r(z) is an operator that preserves 
£%' b (z), but does not preserve V n {z)- By Lemma B.121 p(b) ^ 0. By Proposition ^. 81 
T preserves £ a ' b {z). Hence, Proposition 14 . 51 applies . and we can define polynomials 
qo(z),r (z) so that equations (|5*2|l - (|5*5|l hold. Let ki be defined by equation (J72J), and 
let A(z) denote the difference of the expressions shown in © and (|70|l . By Lemma 
{52) , A(z) is a polynomial such that equation (f?3|) holds. Let T = p(w)D ww + 
q{w)D w + r(w) denote the inversion of T. By covariance and by Proposition 14.41 
T preserves £%' b (w) where 

w = —l/z, a = n(b + ab), b =—1/6. 

Applying the above argument, mutatis mutandi, to T we conclude that A(w) is a 
polynomial such that equation Ij73(l holds with p(z), b, Zi replaced by p(w), b, Wi — 
—1/zi, respectively. By Lemma 14.141 A(w) and A{—l/w) differ by a constant. 
Since both A(w) and A{z) are polynomials, a fortiori, they must both be constants. 
Hence, the potential v(z) has the form shown in l|70() . Furthermore, since A'(b) = 0, 
equation 173() implies 169JI . Finally, since both qo(z) and qo(w) are polynomials and 
since Q72JI holds, Lemma f4. 151 gives equations H68(l and Q81JI. Evaluating the latter 
at z = b, and using (|54|l gives equation l|71(l . □ 



5. Classification of Xi potentials. 

The characterization of Xi potentials achieved by Theorem 14.111 allows to per- 
form a complete classification of these potentials, by using a similar approach to 
the classification of Lie-algebraic potentials. Since the equivalence problem admits 
SL2 covariance and the leading order component of the operator T is a quartic, the 
classification proceeds by considering the 6 possible root configurations of a quartic 
polynomial. As above, we name the classes type I, II, D, Z, III, and N. 

To keep things manageable we focus on non-singular potentials whose domain, 
in the physical variable x, is the entire real line. The type Z change to physical 
coordinate is z = sinhai, so that both the physical and the algebraic domain are the 
entire real line. However, an Xi potential possesses a term of the form l/(z— b), 
where b 6 R, and therefore Xi potentials of type Z are necessarily singular. A 
similar line of reasoning holds for type N potentials. That leaves potentials of type 
I, II, D, and III, which can be obtained by specializing the values of the roots pi in 
Theorem l4.11l In the case of multiple roots, the potential form and the constraints 
are obtained as the appropriate limit of the generic form H70II . 

Some of the resulting potentials also possess multiple algebraic sectors, the anal- 
ysis being essentially the same as for the Lie-algebraic case. 
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5.1. Type I potentials. The Xi potentials of type I are deformations of the gen- 
eralized Lame potentials. We take 

p(z) — Az{z — l)(mz — to + 1), to G (0, 1), 

and denote the roots of p(z) by 

z x = l—— , z 2 =0, z 3 = 1, z 4 = oo. 

TO 

The physical coordinate x is related to the algebraic coordinate z by z — cn(a;|TO) 2 , 
where cn(a;|m) is the Jacobi elliptic function of modulus to. In this way, the physical 
domain, — oo < x < oo, maps onto the interval < z < 1. We now take the 
potential form shown in Theorem 14.111 perform an inversion z = — l/w and to 
ensure that the potential is non-singular demand that k% = k 4 = ±1/4, and that 
b g" [0, 1]. The resulting potential has the following form: 

u{x) = toci( Ci + 1) sn(x\m) 2 +mc 2 (c 2 + 1) - ) ' + (85) 

cm{x\m) z 

P'(P) 2p(6) „ = 1 -2/c- 



cn(x\m) 2 — b (cn{x\m) 2 — b) 2 1 2 
Restating (|68ll l|69|l . the constraints on k\, k 2 are given by 

k\ + k 2 = -n - (k 3 + fc 4 ), k 3 + k A € {-5,0, \} 
,2 m -! ,2 m(b-l) 2 + 2b-l . . 

(oto — m + \y 16 [p — \y b z to 

Let K(m), K'(m) = iK(l — to) denote the real and imaginary quarter-periods of 
cn.(x\m). For b > 1, the equation cn(x|r7i)2 = \fb has a unique, up to addition of 
an integer multiple of 2K' ', solution of the form a; = ip, p € M. A calculation shows 
that the residue of u(x) vanishes at i = ±ip. Hence, up to a constant term, the 
potential can be expressed as 

u(x) = m{ci(ci + 1) sn(a;|m) 2 + c 2 (c 2 + 1) sn(x + (to) |to) 2 (88) 

+ 2 sn(.T + ip\m) 2 + 2 sn(x — ip\m) 2 }. 

For integer values of C\ , c 2 the resulting expression is a finite-gap potential (also 
known as a Picard potential [5]) that was first described in reference [20]. It is 
the simplest known example of a finite gap potential that falls outside the elliptic 
Treibich-Verdier class (the type I Lie algebraic potentials). For more on this see 
reference [17]. 

The algebraic form of the operator is 

„(3)m vW(b) 
T = iL -r 1 J 2 + ^^.h+p'i^Ji+pib)^ + (k 4 + 4TOn) J 6 , (89) 
2 

and it preserves the Xi subspace £° ,fc (z), where 

mk 2 k 3 k 4 



(90) 



mb — to + 1 b 6—1 
The gauge factor is 

p(x) = (cn(ir|TO) 2 — b)~ x dn(x|TO) C2 cn(x|m) C3 sn(a:|TO) C4 , (91) 

where C3, C4 € {0, 1} according to the signs of £3, k 4 . 
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5.2. Type II potentials. The double root is at infinity. The other two roots are 
at zz = and za = 1. Thus, 

p(z) = 4z(l — z), z — — sinh 2 (a: ) 

The non-singularity condition is k 3 = ±1/4 and 6 £ [—1,0]. The potential form 
and the gauge factor are given by 

u(x) = 4£i 2 sinh 4 x + 4^i (£i + 2£ 2 ) sinh 2 x 

+ c 4 (c 4 - 1) scch 2 x+ 8 b f ^ - 4( - 2b 2 V 
(smn x + 6) z Sinn x + 6 

//(a;) =cxp(-^isinh 2 x)(sinh 2 a; + fe) _1 (cosha;) C4 (sinhx) C3 , (93) 

where, as before, Cj = 2fc^ + 1/2 and c 3 e {0,1}. The potential parameters are 
constrained as follows: 

k± + £ 2 = -n-k 3 , fo = ±\, (94) 

t 1 ((2b- 1)4-2^-^^-^=0 (95) 

The algebraic operator 

T= -4J 3 -4(2 6- 1)J 4 - 4 6(6- 1)J 5 - 8^J 6 (96) 
preserves £^ ,fc (z), where 

5.3. Type D potentials. One double root is at infinity; the other double root is 
at zero. Thus, 

p(z) = —z 2 , z — exp(x). (98) 
The potential and gauge factor are 

26 26 2 

u(x) = (he* - 2£ 2 f + (he- + 2£,f + + {&x , (99) 

p(x) = (e x - 6)- 1 exp(£ ie x - £ 3 e - x + [2£ 4 + \) x), (100) 

where 6 < to avoid singularities. The potential parameters are constrained as 
follows: 

2t 2 + 2^ = -n, (101) 
2£ 3 (£ 3 + 2b£ 4 ) - 6 3 £!(26^ - A£ 2 ) = 0. (102) 

The algebraic form of the operator is 

T = - J 3 - 26J 4 - 6 2 J 5 - 2£i J 6 , (103) 

and it preserves the X 4 subspace £^ b (z), where 

a = 4 + | + ^. (104) 
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5.4. Type III potentials. There is a triple root at infinity and a single root at 
Z4 = 0. Thus, 

p(z) = —4z, z = x 2 . (105) 
The non-singularity condition is = ±1/4. The potential and gauge factor are 

u{x) = t 3 2 x 6 - 2£ 2 £ 3 x 4 + A{£ 2 2 + h£ 3 ) x 2 + ^ + (106) 

[x~ — b) x — b 

fx(x) = {x 2 -6)- 1 exp (-^-x i +£ 2 x 2 ^J x Ci , (107) 

where, as before, b < 0, C4 = 2fc 4 + 1/2 and C4 € {0, 1}. The potential parameters 
are constrained as follows: 

£ 1 = -n-k i , h = ±^, (108) 

lb 2 £ 3 2 ~2b£ 2 £ 3 + £ 1 £ 3 +£ 2 2 = (109) 
The algebraic form of the operator is 

T = -4J 4 - 46J 5 + 84 J 6 , (HO) 
and it preserves the Xi subspace £ ®' (z), where 

a = £ 2 -b£ 3 + 1 j. (Ill) 

6. AN EXAMPLE OF AN Xi ELLIPTIC POTENTIAL. 

Examples of Xi hyperbolic and polynomial potentials and their square-integrable 
eigenfunctions have already been presented in [8,9], and also in [11,12] in the context 
of supersymmetric quantum mechanics. Here we present an example of a periodic 
Xi that can be regarded as a modifcation of the well-known Lame potential. The 
form of the potential is 

u{x) = m £(£+ 1) sn(x|m) 2 + (112) 
12to6 2 + 8(1 - 2m)6 + 4m - 4 8b(b - l)(mb - m + 1) 
cn(a;|m) 2 — b (cn(x|m) 2 — b) 2 1 

where m 6 (0, 1), b > 2£/(2£ — 1) and where £, m, b are related by 

/ „ n9 m(b-l) 2 + 2b-l m-l . 

(2£-l) 2 = \ - h (113) 

[b — 1) b z m [bm — m + 1) 

If £ is a positive integer, we set 

fci = - - -£, k 2 = ±- 
4 2 ' 2 4 

so that l|112() becomes a specialization of (|85|l and (|113fl becomes a specialization 
of 18 7|) . A plot of the potential for £ = 4 and various values of m can be seen in 
Figure Q 

From 19 H . the eigenfunctions have the following form 

ij)(x) — (cn(x|m) 2 — b)^ 1 dn(ir|m) C2 cn(x|m) C3 sn(a;|m) C4 p(cn(a;|m) 2 ) 



■2C, 



DAVID GOMEZ-ULLATE, NIKY KAMRAN, AND ROBERT MILSON 





Fig. la m = 0.01 



4 6 8 10 12 

Fig. lb m = 0.33 




Fig. Id m = 0.999 



Figure 1 . Different forms of potential (|112|l for £ — 4 and various 
values of m. 

where Ci — 2ki + 1/2 and p(z) <E £^ ,b \z) with a given by H9()|l . The potential has 
four algebraic sectors. Depending on whether £ is even or odd, the four possible 
algebraizations for each choice of £ are shown in Table El 





£ 


even 








<?odd 






Sector 


n 


C2 


C3 


c 4 




n c 2 


C3 


c 4 


1 


1/2 


1 








1/2 


+ 1/2 








2 


1/2 





1 





1/2 


- 1/2 1 


1 





3 


1/2 








1 


1/2 


- 1/2 1 





1 


4 


1/2- 


1 1 


1 


1 


1/2 


- 1/2 


1 


1 



Table 6. The four algebraic sectors of the Xi elliptic potential 1112|) . 



Having fixed £, the potential has only one remaining free parameter, which we 
take to be m and use (|113f) to determine b, restricted to the condition b > 2£/(2£—l). 
The free parameter m ranges in the interval (0, 1). The two limiting cases are 

(1) in — ► 0, which implies that b — > oo. The potential l|112|l tends to a sin 2 (cc) 
trigonometric potential (see Figure^). 

(2) m — * 1, so that b approaches 2£/(2£— 1). The potential in this limit assumes 
an interesting form of very deep, short range wells on a uniform background 
(see Figure nji). 

For integer values of £, potential (|112|l has exactly £—1 gaps in its band spectrum 
The 2£ — 1 eigenvalues that describe the edges of the allowed and forbidden energy 
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bands (regions of stability or instability) can be calculated algebraically. Note that 
Lame potential m£(£+l) sn 2 (x\) has £+1 gaps in its band spectrum, and the 2£+l 
eigenfunctions corresponding to the band edges are the Lame polynomials. The Xi 
elliptic finite gap potential (|112f) has one gap less than Lame, corresponding to the 
fact that the invariant polynomial space £n' b \z) is a codimension one subspace of 
V n - For instance, when £ = 4, the first three algebraic sectors are two-dimensional 
while the fourth sector is one-dimensional (see TablelHJl. For each of the four sectors, 
the actual eigenvalues and eigenfunctions are calculated by diagonalizing the matrix 
corresponding to the algebraic action shown in (|89|l . As an illustration, for £ = 4 
and to = 0.999 the seven algebraic eigenfunctions corresponding to the band edges 
have been computed and displayed in Figure |3 A look at the energies shows that 
the allowed bands are very narrow, which is natural since the value of m is very 
close to unity. In the limit to — > 1 the period of the potential diverges and the 
bands collapse into pure a point spectrum. 

7. Summary and Outlook 

In this paper we introduced the notion of an exceptional polynomial subspace 
of V n and identified those subspaces as the only ones that can give rise to new 
quasi-exactly solvable potentials in one dimension. For exceptional subspaces of 
co-dimension one (Xi subspaces), we claimed ('Theorem 14.1(1 that £ n is the only 
Xi subspace up to an SL2 transformation. Next, we used the equivalence between 
an arbitrary second order differential operator and a Schrodinger operator to con- 
struct, characterize and classify Xi potentials, which are not Lie-algebraic. The 
characterization is done at the level of the potential invariant, and it involves an 
additional quadratic condition on the residues ki of the quotient q(z)/p(z) of the 
coefficients of the second order operator. As in the Lie-algebraic case, the leading 
order component p(z) of an Xi operator is a quartic and we have used SL2 covari- 
ance to classify them into canonical forms, leading to new quasi-exactly solvable 
families of elliptic, hyperbolic, trigonometric and rational potentials on the line. 
The new characterization of the QES condition at the level of parameters ki in the 
potential invariant allows to analyze the potentials that admit multiple algebraic 
sectors, corresponding to a residual symmetry in the choice of ki such that the 
potential remains unchanged. We have provided a classification of all such cases. 

Future work involves extending this classification into two possible directions: 

(1) Exceptional subspaces of codimension two or higher. Some examples of X2 
subspaces are known to exist, but a full classification is not yet available. 
The equivalence problem is considerably harder than in the Xi case. 

(2) Exceptional subspaces of polynomials in two or more variables. A full clas- 
sification of exceptional Schrodinger in more than one variable operators 
seems unfeasible due to the lack of a constructive solution to the equiva- 
lence problem. Nevertheless, some examples of many-body quasi-exactly 
solvable models are known [2,3], and we do not exclude that new many- 
body problems exist with invariant spaces of exceptional type. 
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FIGURE 2. The periodic potential u{x) in i|112|) for i = 4 and 
to = 0.999 has four gaps in its band spectrum. The seven alge- 
braic eigenfunctions that correspond to the band edges are plotted 
together with their energies and the algebraic sector to which they 
belong. 



